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Abstract 

Motivated by four-dimensional superstring models, we consider the possibility of 
treating the Yukawa couplings of the Minimal Supersymmetric Standard Model 
(MSSM) as dynamical variables of the effective theory at the electroweak scale. 
Assuming bottom-tau unification, we concentrate on the top and bottom Yukawa 
couplings, and find that minimizing the effective potential drives them close to an 
effective infrared fixed line. Requiring an acceptable bottom-top mass ratio leads 
in principle to an additional constraint on the MSSM parameter space. As a by- 
product, we give new approximate analytical solutions of the renormalization group 
equations for the MSSM parameters. 
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1. In a recent paper |l|], the possibility was discussed of treating the Yukawa couplings 
of the Minimal Supersymmetric Standard Model (MSSM) not as numerical parameters but 
as dynamical variables (for similar suggestions, see also This possibility naturally 

arises when one embeds the MSSM into a more fundamental theory, such as supergravity 
or superstrings, where parameters are replaced by vacuum expectation values of some 
singlet scalar fields (moduli), corresponding to approximately flat directions of the effective 
potential. In U, the discussion was mainly restricted to the top-quark Yukawa coupling, 
and it was found that, if the scale Msusy of the explicit MSSM mass terms is of the 
order of the electroweak scale, and if supersymmetry breaking does not induce mass terms 
larger than 0(MsusY 2 /^pianck) in the relevant moduli directions, then minimizing the 
vacuum energy attracts the top-quark Yukawa coupling close to its effective infrared fixed 
point, which is compatible with a top-quark mass in the experimentally allowed range. In 
U, it was assumed that the scale M SU sy, proportional to the gravitino mass m 3 / 2 , also 
corresponds to an approximately flat direction of the fundamental theory, as suggested 
by a certain class of supergravity models ||||. However, we shall see that the result 
on the top-quark Yukawa coupling remains valid even if one assumes (as is often done in 
phenomenological analyses) that Msusy is fixed by some physics at very high scales, which 
allows Msusy to be treated as an input parameter, independent of the Yukawa couplings 
apart from the usual renormalization effects, in the effective field theory at the electroweak 
scale. 

In the present paper, we examine the possibility of dynamically explaining the third- 
generation fermion masses: m t , rrib, m T . To do so, we generalize the considerations of |]J 
to the case where all the third-generation Yukawa couplings are included, still neglecting 
the two light generations. Assuming for simplicity the unification relation hb(Mjj) = 
h T (Mjj), which in first approximation fits the experimental value of the mb/m T ratio H[7j, 
and treating the explicit MSSM mass terms as numerical parameters, we find that in 
this more complicated case the Yukawa couplings are dynamically attracted close to an 
effective infrared fixed line, F(h t , hb) = 0. We also find that minimization of the effective 
potential with respect to the residual variable 8, which parametrizes the infrared fixed line, 
dynamically fixes the ratio hb/ht, and may allow for acceptable values of the mb/m t ratio 
within the residual MSSM parameter space. Our approach leads to the elimination of two 
of the free parameters of the MSSM, ht and hb, in addition to the parameter h T , removed 
by the unification relation hb{Mjj) = h T (Mu). Furthermore, the infrared behaviour of the 
renormalization group equations (RGEs) for the mass parameters of the MSSM is such 
that the latter are severely constrained, with some combinations being driven close to 
effective infrared fixed values. 

The structure of our paper is as follows. In section 2, we review the infrared behaviour 
of the running Yukawa couplings and masses in the MSSM presenting some ap- 

proximate analytical solutions of the RGE for the top and bottom Yukawa couplings^ as 

1 Similar results have been simultaneously and independently obtained in W%. 
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well as for the MSSM mass terms, when both hb and h t are non-negligible. In section 3, we 
present the theoretical motivations that lead us to minimize the effective potential not only 
with respect to the Higgs fields, but also with respect to the top and bottom Yukawa cou- 
plings, taken as independent dynamical variables. In section 4 we describe in some detail 
the results of such a minimization, both analytically and numerically, working for given 
numerical boundary conditions on the MSSM mass parameters at the unification scale 
M\j. We begin by showing that the leading dependence of the effective MSSM potential 
on h t and hb forces the latter, upon minimization, to lie close to the effective infrared fixed 
line. We then proceed to the more subtle problem of minimizing the low-energy effective 
potential along the infrared fixed line. In section 5 we summarize our conclusions, after 
commenting on the case in which Msusy is also taken as a dynamical variable and on the 
case in which the top and bottom Yukawa are constrained dynamical variables. 



2. Neglecting as announced the first two generations, and working as usual in a 
mass-independent renormalization scheme, the one-loop RGEs for the Yukawa couplings 
read |15| 
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where a tAr = h 2 tAT /(4ix), t = \og(M*/Q 2 ) and M v ~ 2 x 10 16 GeV. Two-loop RGEs are 
available |TB|, but we do not need them for our present purpose. 

If in eqs. (|IJ) we neglect a& and a T with respect to a t , which is a good first approximation 
when tan/3 = v 2 /vi <C m t /mb, and take Q ~ Msusy ~ fnz, then the effective infrared 
fixed point for the top-quark Yukawa coupling is ||[17| oc t — (8/9) «3, corresponding to a 
running top-quark mass m t ~ (195 GeV) sin (3, from which the pole mass can be extracted 
by including the standard QCD corrections, mP° le = m t (m 4 ) [l + 4a 3 (m t ) / (37r) + . . .]. When 
ab and a T are not neglected, which is the case of interest for the present paper, eqs. (P 
have a more complicated infrared behaviour. Choosing for simplicity a 



a; 



in order to 

have a more manageable two- variable problem, the resulting infrared structure is displayed 
in fig. 1. We choose random boundary conditions satisfying the constraint 



2au < a¥ + a% < 1 



(2) 



with au — 1/25, corresponding to the dots in the region of the (a^,a^) plane shown 
in fig. la. We then solve numerically the RGEs of eq. ([!]) at the representative scale 
Q = 200 GeV, to obtain a t , ab and a T . The resulting region of the (h t , hb) plane is shown 
by the dots in fig. lb. The important aspect to be stressed is the focusing effect due to the 
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infrared structure of the RGE: a relatively wide region of the (h^ , h^) plane is mapped 
into a very thin region of the (ht, hb) plane, clustering around an 'effective infrared fixed 
line'. Another effect, clearly visible in fig. 1, is the existence of some special points along 
the effective infrared fixed line. If we look at the density of points in the (h t , hb) plane, 
corresponding to a uniform distribution in the (h^ , h^) plane, we can clearly see that the 
point h t = hb is an attractor, whereas h t = and hb = are repulsors. 

For practical purposes, we now introduce some approximate analytical formulae for 
ht and hb, which can be useful to parametrize the effective infrared fixed line. If in the 
RGEs for a t and we neglect the terms proportional to a T and to a', and we define the 
auxiliary variables 

p = y/h$ + hl, tan0 = ^, (3) 
after some calculations we can write the solution as 

/ 'fori 2 29) /(sin 2 20 v ) 1 J_F 

P 2 pi e + ^e' [) 

^'"'.E, (5) 



p 2 (sin 29) 12 / 5 p 2 / (sin2^^ 12/5 



where 



(COS 2^) 7 /5 ( cos 20^)7/5 

f(x) = 2 T X (1/2, 1, 11/5 ; x) = \ I" ds{\ - sf'^l - sx)- 1 ' 2 , (6) 
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The same result was independently obtained in ref. |^4| . Since the behaviour of the 
hypergeometric function / will play an important role in the considerations of section 4, 
for illustration we plot in fig. 2 /(sin 2 29) and <i/(sin 2 29) / d9, as functions of 9. From eq. (|), 
we can see that the renormalization of the quantity in the first member is multiplicative. 
On the other hand, 1 < /(sin 2 29) < 12/7, and at scales Q ~ 200 GeV it is (still neglecting 
a' effects) E ~ 10 and F ~ 220, so that the second member of eq. (^) is dominated by 
3F/(8ii 2 E) for sufficiently large values of pu (pu ~ 0.5). This defines an effective infrared 
fixed line at any scale Q close to the electroweak scale, parametrized by 

h t = Pir{9) cos 9, h b = Pm{9) sin 9, (8) 

where 

Pm(9)^ ] j^f(sm 2 29). (9) 

To further improve the approximation of eq. we can introduce some constant shifts to 
fit the corrections due to a' and a T effects, and get 

ht = 0.015 + pm(9) cos 9, h b = -0.045 + p JR (8) sin 9 , (10) 
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where of course the range of variation of 9 should be modified accordingly. It is useful to 
compare our formula with the one previously derived in [[13], which, after correcting for a' 
and a T effects as in eq. ([TUP, reads 



(h t - 0.015) 12 + (h b + 0.045) 12 = [j^J . (11) 

Figure lb compares the exact numerical solutions of eqs. (|1|), represented by the dots, with 
the approximate analytical solutions of eqs. ( PH| ) and (|TT|), represented by the solid and 
by the dashed line, respectively. We can see that both formulae are good approximations 
for ht <C h t or h t <^ hb, whereas eq. ( |10D is a better approximation for hb ~ h t . 

An essential ingredient in the study of the MSSM effective potential is the solutions to 



the RGE for the MSSM mass parameters [18]. Exact analytical solutions of the one-loop 



RGE are known [19] in the case of negligible hb and h T . Some approximate analytical 



solutions have also been obtained recently for the special case h b = h t |12| . We have 
improved the existing formulae by constructing approximate analytical solutions valid for 
any value of h t and hb-, and including the most important a' effects, but still neglecting 
a T effects. Their explicit form is given in the Appendix. From our formulae one can 
easily rederive some known relations valid at special points on the effective infrared fixed 
curve: A t /m 1/2 ~ H 8 - H A /2 ~ 1.5 and Am 2 ~ (3/2)m 2 + (l/2)(H 2 - Hj/A)m 2 1/2 ~ 
(3/2)rng + (l/2)(6.3)m 2 /2 for B = Q, A t /m 1/2 ~ A b /m 1/2 ~ 1.5 and Am 2 ~ Am 2 ~ 
(9/7)m 2 + (3/7)(6.3)m 2 /2 for 6 = tt/4 0. 



3. We now present the theoretical motivations that lead us to minimizing the MSSM 
effective potential not only with respect to the Higgs fields, but also with respect to the 
top and bottom Yukawa couplings, taken as independent dynamical variables. 

In a generic supergravity model, masses and couplings are field-dependent functions. 
This should be kept in mind when considering both the dimensionless and the dimen- 
sionful parameters of the MSSM, seen as the low-energy effective theory of an underlying 
supergravity model. For each given MSSM parameter, if the scalar fields that control it 
are frozen to their VEVs by sufficiently heavy mass terms, then such a parameter can be 
treated as constant, apart from standard renormalization effects, when discussing the dy- 
namics at the electroweak scale. If, however, after integrating out the superheavy degrees 
of freedom, some extra singlet scalar fields are left, with no renormalizable couplings to the 
MSSM fields and masses C(MsusY 2 /^fpianck) or smaller, then MSSM quantum corrections 
can play a role in the determination of their VEVs, and the corresponding MSSM param- 
eters should be treated as dynamical variables of the effective theory at the electroweak 
scale. 

Consider for example the class of supergravity models whose field content splits into 
an 'observable sector', containing the MSSM states (and possibly others) and a 'hidden 
sector', coupled to the observable sector only via interactions of gravitational strength. 
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The situation of interest to us can be realized in the special subclass of models [|,|J that 
exhibit, in their hidden sector, some approximately flat directions of the classical poten- 
tial, associated to some 'moduli' fields. Such degeneracy of the classical vacuum is in 
general removed by quantum corrections, including the perturbative ones if supersym- 
metry is spontaneously broken. We would like to envisage here the possibility that the 
potential along some of these flat directions does not get large quantum corrections from 
the superheavy sectors of the theory. Then we need to minimize the effective potential at 
the electroweak scale to fix some moduli VEVs and to determine those MSSM low-energy 
parameters which carry a non-trivial dependence on such moduli. 

The above possibility is supported by the general structure of four-dimensional su- 
perstring models ||20|| , where all the low-energy parameters, in particular the Yukawa 
couplings, are dynamical variables depending on some moduli VEVs. Indeed, some in- 
teresting superstring solutions could give rise, in the low-energy limit, to spontaneously 
broken N = 1 effective supergravities of the type considered above. At the classical level, 
gauge and Yukawa couplings are related |21| by a string super- unification condition: 

^ = — . (12) 

In eq. (|T2"D, a str is the coupling constant associated with the string loop expansion: already 
at the classical level, this is not a numerical parameter but a dynamical variable, related 
to the VEV of the dilaton field (S + S) by 2ira str = (S + S)" 1 . For string solutions with 
unbroken supersymmetry, a str is a flat direction, not only classically but at all orders in 
the string perturbative expansion. In the case of the gauge couplings (i = 1, 2, 3 for the 
factors of the standard model gauge group), the coefficients ki are constants depending on 
the particular string solution. In the following, we shall have in mind the class of string 
solutions for which, with the usual normalization convention gi = JhJZg' ', it is k 3 = k 2 = 
k\ — 1. We shall also assume that some non-perturbative effects break spontaneously N = 
1 supersymmetry and fix the VEV of the dilaton field, in such a way that gauge couplings 
are not dynamical variables at the electroweak scale. In the case of the Yukawa couplings 
(i = t, b, r for the third-generation ones to be considered here), at the string classical 
level the coefficients ki typically are exponentially suppressed or of order unity, as can be 
easily checked in many explicit examples. For instance, in free fermionic constructions the 
non-vanishing tree-level Yukawa couplings correspond to moduli-independent coefficients 
ki = 2. In other classes of string solutions one can have moduli-dependent ki coefficients 
for the tree-level Yukawa couplings. 

Even when the tree-level moduli-dependence of the Yukawa couplings is identical to 
that of the gauge couplings, string-loop corrections P2| to the low-energy effective action 
can introduce additional, non-universal moduli dependences. One can then envisage, as 
already discussed in |TJ, various possible situations. 

A first possibility is that, after the inclusion of string- loop corrections and of possible 
non-perturbative effects associated with supersymmetry breaking and with the stabiliz- 
itaion of the dilaton VEV, the Yukawa couplings have no residual moduli dependence. In 
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this case, the non-vanishing ones will still obey a superunification condition of the form 
(|T2|), with ki = 0(1). In particular, the top and bottom Yukawa couplings will fall in the 
domain of attraction of the infrared fixed curve discussed in section 2. 

A second possibility is that for some Yukawa couplings there is a residual moduli 
dependence along some approximately flat directions. In the following section, we shall 
assume that this moduli dependence preserves the unification relation h b = h% , but 
allows the top and bottom Yukawa couplings to be treated as independent variables of the 
effective theory at the electroweak scale. Of course, more complicated situations could also 
arise, for example that there be only one independent flat direction in moduli space. In 
this case, still assuming for simplicity h% = h 1 ^ , the allowed range of variation for h% and 

would be restricted to a certain curve of the (h^ , h%) plane, and this should be taken 
into account when minimizing the low-energy effective potential. We shall temporarily 
disregard this last possibility in the following section, but we shall come back to it in the 
concluding one. 



4. If, as suggested by the models of ref. ||, there are no quantum corrections to 
the vacuum energy carrying positive powers of superheavy scales, the one-loop effective 
potential of the MSSM can be written as V = Vq + AV, where 



Vq = m\v\ + m\v\ + 2m\v \V 2 + 9 — 9 (vj - vl) 2 + r] , (13) 



and 



mf 



AV =^ E(-l) 2Jl+1 (2^ + l)mf (log ^ - -) ■ (14) 

In eq. (0), V\ and v 2 are the neutral Higgs vacuum expectation values, and ml = m 2 Hi +fi 2 , 
m\ = m 2 H2 + fi 2 , m 2 = Bfi are mass parameters. The 'cosmological term' rj = fjm\i 2 takes 
into account, as in |1] but in a slightly different notation, the contributions to the low- 
energy effective potential that do not depend on the MSSM fields. Given a set of boundary 
conditions at Mjj, the mass parameters m 2 ,m 2 ,m 2 have an implicit dependence on the 
top and bottom Yukawa couplings via their RGEs, as illustrated by the approximate 
analytical solutions given in the Appendix. A similar implicit dependence is present for 



the parameter rj, whose renormalization-group evolution was studied in Q. In eq. (|i4|), 
rrii and Jj are the tree-level field-dependent mass and the spin for each particle % in the 
MSSM spectrum. Notice that AV has an explicit dependence on h t and hb only via the 
top, bottom, stop and sbottom squared masses, 

m t = h 2 t vl , m 2 b = h 2 h v\ , (15) 
m h 2 = h t v 2 + + a ( Vl ~ v<2 > 
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As announced, we shall discuss here the minimization of the MSSM one-loop effective 
potential not only with respect to the Higgs fields, V\ and v 2 , but also with respect to the 
top and bottom Yukawa couplings, ht and hj,, treated as independent dynamical variables. 
We would like to stress once more the importance of the cosmological term rj in the RG- 
improved tree-level potential of eq. (|I3|). This term is usually neglected because it does not 
depend on v\ and v 2 \ hence it does not play any significant role in the minimization with 
respect to the Higgs fields. In our case, however, this term must be included, since, given 
a boundary value r] = T](Mu), tj has an implicit dependence on the Yukawa couplings via 
its renormalization group evolution: neglecting 77 would create an artificial dependence of 
the effective potential on the renormalization scale Q. 

With the above comments in mind, we can proceed to the minimization of the one- 
loop effective potential with respect to the Higgs fields and the top and bottom Yukawa 
couplings, given a set of boundary conditions (mo, mi/2, A , B , p ; 770) ■ For convenience, 
we work with the polar coordinates p and 9 already introduced in section 2, and we proceed 
in two separate steps. First, we fix 9 to an arbitrary value, and we minimize the potential 
with respect to vx, v 2 and p. We find that, for any given value of 9, the value of V at its 
minimum with respect to v\ and v 2 gets smaller and smaller as p increases, until p reaches 
its maximum allowed value, pm(0), corresponding to a point on the effective infrared 
fixed line. This result has been tested numerically for many different values of 9 and of 
the boundary conditions on the free parameters, using the full one-loop effective potential 
of eqs. ( |13D and (14) . We also verified that in most cases, with an appropriate choice of 
the renormalization scale, Q 2 ~ mj , the minimization with respect to h t and hj, is 
dominated by the V contribution: this extends the results obtained in |^3| for the usual 
minimization with respect to v\ and v 2 . The mechanism of attraction towards the effective 
infrared fixed line can be understood semi-analytically in sufficiently simple cases, as we 
shall now discuss on an example. 

Consider the toy version of the MSSM corresponding to mo = Aq = Bq = po = 0, 
vx = 0, hb = 0, but with mx/2 and ?7o both taken as fixed numerical inputs and not as 
dynamical variables. In this case, eq. (Op simplifies to 



Vq = mlv 2 + 



g 2 + g' 2 



v 2 + r). 



(18) 
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Assuming m\ < 0, as needed for SU(2) x U(l) breaking, and minimizing with respect to 
v 2 , we find Vq]^^ = -2m\j(g 2 + g' 2 ) + 77, and finally 



dx 



1 -2m\ dm l ,dv nq x 



, v 2n a 2 + a' dx dx 

V2 = (v 2 ) Z 



where x = a t /aj R = 2nEa t / (3F). The Yukawa-coupling dependence of m 2 , can be easily 
understood by specializing the general formulae given in the Appendix, 

ml = [C + Ax + Bx 2 )m\ /2 , (20) 

where, at scales Q of the order of the electroweak scale and in the notation of the Appendix, 
C = Ci/4 + C 2 ~ 0.5, A = -# 2 /2 ~ -5, 5 = #f/8 ~ 2. We can check the well-known 
fact that, in the physical region x < 1, we always have Ac+.Ba; 2 < 0, which leads to m\ < 
for sufficiently large values of x. If the 77-dependent part of eq. (|19D can be neglected, we 
can already state that in the case under consideration x is driven to x — 1. However, we 
know from that 77 (slowly) increases for increasing x, thus a quantitative comparison of 
the two terms in eq. flT9| ) is necessary. From the RGE for rj and a t , we obtain 

dv 1 (A + Bx)(D + Ax + ^x 2 ) 4 

dx" 2tt G + J(l-x) mi / 2 ' ( j 

where, at scales Q of the order of the electroweak scale and in the notation of the Appendix, 
D = 2C 3 + 2C 2 + 13Ci/18 ~ 11.3, G = 16a 3 /3 + 3a 2 + 13a'/9 - 6«P ~ 0.01, J = 6«P ~ 
0.57. From this one can easily verify that the 77-dependent part of eq. (|19"D is indeed 
negligible for all values of x < 1 . 

Having obtained the result that, for any given value of 9, minimization with respect to p 
invariably leads to p = pm(9), we can now restrict our attention to top and bottom Yukawa 
couplings constrained along the effective infrared fixed line, and minimize the effective 
potential of the MSSM with respect to the residual angular variable 9 (in addition to the 
usual variables v\ and v 2 ). Numerical investigations show that, depending on the chosen 
boundary conditions for the mass parameters, different structures may appear. A typical 
situation is illustrated, for a representative parameter choice, in fig. 3: it corresponds to a 
trivial minimum for 9 = 0, i.e. to vanishing bottom and tau tree-level masses. As we shall 
discuss later, this type of structure can be rescued by some constraint on the moduli space 
that forbids the boundary condition 9jj = 0. In this case, the low-energy 9 just relaxes 
to its minimum allowed value, 9 m i n 7^ 0, and a hierarchy < mb/nh "C m t can emerge 
even for values of tan/? close to 1. In our numerical investigations, we were not able to 
find unconstrained non-trivial minima for 9 ^ 0, corresponding to universal boundary 
conditions on the mass parameters and a particle spectrum compatible with the present 
experimental data. Establishing whether realistic solutions of this kind can be obtained 
or not would require further investigations. If such solutions do exist, minima close to 
9 = it/A would be favoured by the peculiar behaviour of the function f (sin 2 29) and by 
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the focusing effect of the RGEs. When 9 ~ 7r/4, acceptable values for m t and mj can 
be obtained only for tan/3 ~ mt/m&. This situation can be realized either by selecting a 
strongly restricted region of parameter space or by allowing some violation of universality 
in the boundary conditions for the mass parameters. 



5. Motivated by four-dimensional superstring models and their effective supergravity 
theories, we examined the possibility of treating the Yukawa couplings of the MSSM as 
dynamical variables at the electroweak scale. In particular, we concentrated on the Yukawa 
couplings of the third generation, neglecting the two light generations and assuming the 
unification relation h% = h% ■ We have found that, treating (h^h^) as independent 
variables, minimization of the one-loop MSSM effective potential attracts (h t , h^) to an 
effective infrared fixed line. This general feature allows the elimination of one of the free 
parameters of the MSSM, and leads to the generic prediction 

8 32 

-a 3 <a t + a b < —a 3 , (22) 

which can be further improved by including a T , a', higher-loop and threshold corrections. 
In terms of the top and bottom quark running masses, the above prediction reads 

{Ml R f <^% + ^%< ^(M«) a , (23) 
sin p cos^ p 7 



where M/ B ~ (4/3)^/a 3 /(a 2 + o>')m z — 195 GeV. This result can be translated into a 
relation involving the pole top and bottom masses by straightforward inclusion of some 
finite MSSM one-loop corrections, dominated by standard QCD effects. The ratio otb/ott 
is also determined by minimization, but its actual value at the minimum depends on the 
free mass parameters of the MSSM (m , m^, A , B , /x ; rj ). The number of the MSSM 
free parameters is further reduced by one, but no generic prediction can be made in the 
absence of a theory of the mass parameters. 

Our results were obtained under two important assumptions. First, the overall scale 
^susy of the MSSM mass parameters, proportional to the gravitino mass, was not taken 
as a dynamical variable but as a fixed numerical input. Second, the two Yukawa couplings 

and h% = h 1 ^ were considered as independent variables in the minimization. We would 
like to conclude our paper by commenting on the effects of relaxing each of these two 
assumptions. 

When, as in ]1J, M SU sy is also considered as a dynamical variable, one obtains an addi- 
tional constraint on the MSSM mass parameters, coming from the minimization condition 
with respect to the gravitino mass, which sets the overall MSSM mass scale, 



., dVi StrJW 4 

1 + ^ = - (24) 
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It was shown in QXJ] that the above equation allows for the dynamical generation of the 
desired Msusy /^Planck hierarchy in a large region of the parameter space. The main quan- 
titative result on the dynamical determination of the Yukawa couplings, i.e. the generic 
attraction of a t + towards the effective infrared fixed line, remains the same. However, 
the determination of the ratio a&/at as a function of the boundary conditions for the 
residual free mass parameters would require a separate analysis, since the minimization 
condition with respect to 7713/2, eq. (H), induces a non-trivial dependence of the overall 
mass scale of the effective potential on the angular variable 9 parametrizing the effec- 
tive infrared fixed line: this might allow for an easier generation of phenomenologically 
acceptable non-trivial minima at 9 ^ 0. 

Another possibility is that the top and bottom Yukawa couplings at the unification 
scale, h% and h}( , are not independent but constrained by some functional relation, corre- 
sponding to a curve in the {h^ , h%) plane, such as those shown in fig. lc. Such a possibility 
can occur if the moduli dependences of h!( and are correlated, and correspond to a 
single independent flat direction in moduli space: in this case also the number of inde- 
pendent minimization conditions has to be restricted accordingly. However, the generic 
phenomenon of attraction towards the effective infrared fixed line will persist also in this 
case, as long as the constraint on (h^ , h%) allows for some points with sufficiently large pu 
to fall in its domain of attraction. If the constraint at the unification scale allows for all 
possible values of 9\j within the domain of attraction of the effective infrared fixed line, 
then such a curve is mapped into the entire infrared fixed line, and, as far as low-energy 
Yukawa couplings are concerned, minimization under the constraint gives exactly the same 
result as unconstrained minimization. It might well be, however, that either the range of 
variation of 9jj is restricted, or a sufficently large value of pu is allowed only for certain 
values of 9u, as is the case for the dot-dashed and solid curves in fig. lc, respectively. In 
this case the minima of the low-energy potential will still lie on the infrared fixed line, but 
minimization with respect to 9 must take into account the bounds set by the constraint 
at the unification scale, as apparent from the corresponding curves in fig. Id. The only 
case in which the constrained minimum does not lie along the effective infrared fixed line 
corresponds to a curve in the , h%) plane that does not allow for sufficiently large values 
of pu- we regard this last situation, exemplified by the dashed lines in figs, lc and Id, as 
extremely unlikely, given the fact that the tree-level string Yukawa couplings are typically 
of the order of the unified gauge coupling, which falls already in the domain of attraction 
of the effective infrared fixed line. Several possible constraints on the Yukawa couplings at 
the unification scale were recently conjectured in ||, but without dwelling into a possible 
string origin. We have argued here that the detailed form of these constraints may or may 
not be relevant for the determination of the low-energy Yukawa couplings. The present 
understanding of the moduli space of four- dimensional superstring models, with its gener- 
alized duality symmetries, indeed suggests the possible existence of such a constraint, but 
does not allow to single out a specific form for it in a model-independent way. 
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Appendix 

We present here some approximate analytical solutions to the RGE for the MSSM mass 
parameters, which include h t and hb effects for any values of the latter, but still neglect h T 
effects. We assume for simplicity universal boundary conditions at the unification scale: 

M 3 (M V ) = M 2 (Mu) = M x (Mu) = m 1/2 , (25) 

m 2 Qa (Mu) = m 2 us (Mu) = m^(M a ) = m 2 La {M v ) 

= m% a (Mu) = m 2 Hi (Mu) = m\ 2 (M v ) = m 2 , (26) 

A U (M U ) = A D (M U ) = A E (M U )=A , B(Mu) = B 0} ^{M v ) = ^. (27) 

To parametrize the dependence of the relevant mass parameters on the top and bottom 
Yukawa couplings, we introduce the auxiliary variables 



3F / \ V 3F 



On the infrared curve of eq. (§), we can write x = /(sin 2 29) cos 2 9 and y = /(sin 2 29) sin 2 9. 
To include the most important a' effects, in such a way that our approximate formulae 
are optimized for x £ y, we define 

16 13 ft 

E = Zi Zz Z Z x 99 , F= E(t')dt' , (29) 

Jo 

where (i = 1, 2, 3) 

* - (i + |t) 1 , (30) 

33 

b 3 = -3, b 2 = l, h = —. (31) 

5 

The low-energy mass parameters with a non-trivial dependence on x and y can be written 

as 



and 



, 1 „ „ „ , , 1 

m 



q 3 = m l + (gg^i + °2 + C 3 )m 2 1/2 - - (Am 2 + Am 2 ) , (32) 

(33) 
(34) 
(35) 
(36) 





= ™o + ( 


|d + C3K/2 - 


2 A 2 

- -Am 

3 2 




= m 2 + ( 


-d + C 3 )m 2 /2 - 


2 A 2 

- -Am 

3 1 


m ii 


= m 2 + 


(id + C 2 )m 2 /2 


- Am 2 , 




= m 2 + 


(jd + C 2 )m? /2 


-Am 2 , 
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where 



/ \ 3 / 7 

2 _ 2 / a t a b \ 32/21 7 -3/7 7 -l/231 /o 7 A 

B = B - ^A (x + y) + m 1/2 [H 9 - ^H 4 (x + y)\ , (38) 
A t = A (l -x- l -y)+ m 1/2 [H 8 - ^H 4 (x + ±y)] , (39) 

A b = A (l --x-y) + m 1/2 [H 8 - \h 4 {^x + y)\ , (40) 

1 1 

A T = A (l - -y) + m 1/2 {H 10 - -H A y) , (41) 

d = ^-(1 - 2?) , C 2 = \{l-Zl), C s = -1(1-2%). (42) 

The only two independent quantities entering the solutions for the soft scalar masses are 

3 1 1-1 

Ami = 2 m oy + 2 A ° y ^ ~ a ( x >y)y\ + A o) + ^^[^ - ^ a(x,y)Hjy] , (43) 

3 1 11 

Awij = t^o^ + 2^[1 - a(x, y)x] {H A m l/2 + A ) + -m?/ 2 z[#2 - ^ a(x, y^x] , (44) 

where 

7/[7-x s V] +23 
°(^I/) = ^^ ( 45 ) 

is a suitable interpolating function, such that a(x, 0) = a(0, y) = 1, a(:c, y = x) = 7 /6, and 
ff 2 = |«f 8 , -ff 2 = |iff 8 , H, = 2 (t| - l) , (46) 

4^(^ + 3* + ^), i^^ + fz,). (48) 

The above formulae have been tested numerically by comparing them with the exact 
numerical solutions of the one-loop RGE. If one compares with exact numerical solutions 
neglecting a T effects, our approximate results are correct with less than 3% error. When a T 
effects are included in the comparison, the error of our formulae grows up to a maximum 
of 10%. 
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Figure captions 



Fig.l: Mapping of the {h% ,h%) plane into the (h t ,h b ) plane, for Q = 200 GeV, h% = hF T , 
M\j = 2 x 10 16 GeV, ajj = 1/25. In (b), the dots correspond to the exact numerical 
solutions of the one- loop RGE of eqs. ([!]), for the boundary conditions given in (a); 
the solid line corresponds to the approximate analytical solution of eq. (fCTf), and the 
dashed line to the approximate solution of eq. ([Til) . I n ( c ) and (d), we show how 
some possible constraints, corresponding to curves in the (h^ , h b ) plane, are mapped 
into corresponding curves in the (h t , hb) plane. 

Fig. 2: The function /(sin 2 29) and its derivative df (sin 2 29) / d9 . 

Fig. 3: V{9) for a representative choice of the boundary conditions (m ,mi/2, A , B , /i ) 
and for i] = 0. For convenience, a non- universal contribution S to the boundary 
condition on has been allowed: m? Hi (M\j) = m 2 , + 5 2 . 
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